Inert Doublet Model and LEP II Limits 
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The inert doublet model is a minimal extension of the standard model introducing an additional 
SU(2) doublet with new scalar particles that could be produced at accelerators. While there exists 
no LEP II analysis dedicated for these inert scalars, the absence of a signal within searches for super- 
symmetric neutralinos can be used to constrain the inert doublet model. This translation however 
requires some care because of the different properties of the inert scalars and the neutralinos. We 
investigate what restrictions an existing DELPHI collaboration study of neutralino pair production 
can put on the inert scalars and discuss the result in connection with dark matter. We find that 
although an important part of the inert doublet model parameter space can be excluded by the LEP 
II data, the lightest inert particle still constitutes a valid dark matter candidate. 
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I. INTRODUCTION 

The inert doublet model (IDM) extension of the par- 
ticle standard model has recently attracted attention 
within the particle physics and cosmology community. 
Despite its simplicity - a two Higgs doublet model with 
an imposed unbroken Zi parity - it displays a remarkably 
rich phenomenology. After being touched upon already 
in the 1970s [l| the model has in later years been studied 
within various contexts such as light neutrino generation 
0,1, lepto genesis H, Higgs phenomenology P, Q , im- 
proved naturalness [J, |6( , electroweak symmetry breaking 
Q , and as a dark matter candidate @, H, H, H, Q . 

Besides the standard model-like Higgs boson, h, the 
IDM contains three additional scalar fields: two neutral, 
H° and A , and one charged, H^. Because of the exact 
Z2 parity none of the new particles have direct (Yukawa-) 
couplings to fermions and are hence called inert scalars 
(or sometimes even inert Higgs bosons [3, Q ) . Because 
the new doublet still couples directly to h and the gauge 
bosons (and does not generate standard model masses) 
some have preferred to instead denote the model the dark 
scalar doublet model 

The Z2 parity also ensures the stability of the lightest 
inert scalar, which hence, if it is either H° or A , con- 
stitutes a good dark matter candidate belonging to the 
class of weakly interacting massive particles (WIMPs). 
For masses below 80 GeV it can leave a relic abundance 
in agreement with the WMAP constraints [l(| [H| and 
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hence explain all the dark matter 0, H, Q . 

Another striking property of the IDM is that it can al- 
low for a heavy (~500 GeV) standard model Higgs boson 
h without violating data from electroweak precision tests 
[J] . Although the direct Higgs searches are not in conflict 
with the 95% C.L. upper Higgs mass bound set by the 
precision tests (assuming only standard model physics), 
they are well above the indicated central value [lj] • If up- 
coming accelerator searches fail to detect a light (< 200 
GeV) Higgs boson, the IDM may constitute an attractive 
explanation. 

Although the IDM (in contrast to e.g. supersymmctric 
models) seems to lack any deeper fundamental motiva- 
tions, its strength instead lies in its simplicity, and the 
model can be regarded as an archetype for more com- 
prehensive extensions of the standard model (like [HI or 
[lj|). Just like generic spin- 1/2 and spin-1 WIMPs could 
be represented by the supersymmctric neutralino and the 
first Kaluza-Klein excitation of the photon, respectively, 
the lightest inert particle could be used as a representa- 
tive of spin-0 WIMP dark matter candidates. 

Experimental bounds and signatures of the IDM have 
been investigated in several papers. Studies show that 
the WIMP-nucleon scattering cross section in general 
lies below the sensitivity of current deep underground 
direct detection experiments 0, HI ■ The prospects for in- 
direct detection via gamma-rays produced in dark matter 
self-annihilation pro cesses in the galactic halo look more 
promising [1, I9L Il5j. especially since inert scalar dark 
matter has the ability to produce extraordinarily clear 
spectral lines @ . LEP I data on the width of the Z boson 
force the sum of the H° and ^4° masses to be larger than 
the Z mass [H, Q , while the LEP electroweak precision 
tests put constraints on the inert scalar mass splittings 
as a function of the h mass Q. Moreover, detection at 
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the upcoming LHC seems feasible 

What has been missing so far in the investigation of 
the IDM is a closer study in connection with the existing 
LEP II data. So far the LEP II limits on the IDM are only 
rough estimates based on direct usage of production cross 
section limits from existing analyses of the minimal su- 
persymmetric standard model (MSSM) 0, [H| . Although 
this may be suitable as a first approximation, there are 
certainly some important differences between the IDM 
and the MSSM which threaten to alter the conclusions. 

In this paper we take a closer look at the limits LEP 
II data can set on the IDM. We make use of an existing 
DELPHI collaboration study of neutralino pair produc- 
tion [l6| where detailed signals and backgrounds have 
already been simulated. By generating H°A° events and 
mimicking the cuts performed on the dataset in (l6j , we 
can estimate upper limits on the corresponding produc- 
tion cross sections. These limits are then used to con- 
strain the masses of the inert scalars. 

The outline of the paper is as follows. In Section II 
we present the IDM setup in more detail and discuss its 
theoretical and experimental constraints. In Section III 
we motivate the need for a more careful LEP II analy- 
sis, and introduce the method used here to constrain the 
IDM parameter space. In Section IV we present the re- 
sults, and discuss those in connection with dark matter. 
Section V includes some final remarks and conclusions. 
A short summary of this paper is given in Section VI. Ap- 
pendices A and B provide detailed lists of the cuts used 
for the acoplanar jets and acoplanar leptons selections, 
respectively. 



II. THE INERT DOUBLET MODEL 
A. Description 

Although electroweak symmetry breaking can be 
achieved with a single Higgs doublet, a more complicated 
Higgs sector is not excluded. Since problems like dark 
matter and neutrino mass generation seem to require new 
physics beyond the standard model, minimal extensions 
could be considered as attractive models compatible with 
the principle of Occam's razor. 

A simple and well known extension of the standard 
model Higgs sector is the two Higgs doublet model 
(THDM), in which there exist two scalar doublets, Hi 
and H 2 ■ In order not to be in conflict with the tight con- 
straints on the magnitude of flavor changing neutral cur- 
rents a Z 2 parity, under which Hi —>Hi and H 2 ^ —H2, 
is imposed. The most general renormalizable and CP 
conserving potential of such a model is 

V = vi\Hi\ 2 + n\\H 2 f + A^l 4 + X 2 \H 2 \ 4 



h\Hi\ 2 \H 2 \ 



X A \H{H 2 \ 2 + A 5 i?e {H\H 2 f ,(1) 



of course depends on the values of the parameters in the 
potential [l[ . Even if one usually assumes that both dou- 
blets develop vevs, there is a priori nothing which re- 
quires that phase. 

The IDM is a THDM where only one of the Higgs 
doublets, Hi, acquires a vev, v. Hence, Hi closely cor- 
responds to the ordinary standard model Higgs doublet. 
Moreover, the IDM belongs to the class of type-I THDMs 
in which all standard model fields are taken to be Z 2 
even. As a consequence no Yukawa terms including H 2 
are allowed by the symmetry. The nonexistence of a vev, 
v 2 , for H 2 leaves the imposed Z 2 parity unbroken and 
ensures the absence of mixing between the components 
of Hi and those of H 2 . Thus, the fields belonging to H 2 
are inert in the sense that they do not couple directly to 
fermions, and the lightest of them is automatically sta- 
ble. It should be noted that the IDM is not the t>2 ^ 
limit of a THDM in which both Higgs doublets develop 
vevs. 

After giving masses to the gauge bosons, Hi has one 
physical degree of freedom left: the real scalar field h. 
Since h closely resembles the Higgs particle of the stan- 
dard model it will here be called the standard model 
Higgs boson. In addition, H 2 includes the neutral CP- 
even H°, the neutral CP-odd A , and the charged iJ ± 
inert scalars. The masses of the particles are (at tree 
level) given by 



m 2 H± 



A, 



-2/4 

f4 + (a 3 

(4 + (A 3 + A 4 
\i\ + X 3 v 2 , 



A 5 )« 2 , 
X 5 )v 2 , 



(2) 



where v = m/ l /v / 4Ai is the vev of Hi. (The measured 
gauge boson masses determine v f» 175 GeV.) 

If either of the electrically neutral H° or A° constitutes 
the lightest inert particle, it could be a good WIMP dark 
matter candidate. By assuming m H a < m^o we will from 
now on have H° as our potential dark matter particle, 
although the roles of H° and A are in general inter- 
changeable. 

The mass difference between the ^4° and H° is fre- 
quently appearing in calculations, and for later conve- 
nience we define 



Am = mjifi — mjj° ■ 



B. Constraints 



(3) 



The requirements of vacuum stability and pcrturbativ- 
ity set theoretical constraints on the model. If and only 
if 



Ai, 2 > 
A3, A 3 + A 4 - |A 5 | > -2VA1A2" 



(4) 



where (A 2 and \ are real parameters. Whether the Higgs 
doublets develop vacuum expectation values (vevs) or not 



the potential V of Eq. (Q]) is bounded from below, en- 
suring a stable vacuum. We will follow [J] and adopt 
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A§ + (A 3 + \a) + X 2 5 < 12Xl 
A 2 < 1 



(5) 



as sufficient conditions for a perturbatively well-behaved 
model. 

Observational constraints come from direct and indi- 
rect detection experiments, measurements of the cosmic 
microwave background, and accelerator searches. 

The absence of WIMP-nucleon scattering signals in ex- 
isting direct detection experiments disfavors models very 
close to having a Peccei-Quinn symmetry mno = tuao 
[J]. Other direct detection constraints on the IDM pa- 
rameter space arc not considered in this work [4,8] . 0, Q . 

Indirect detection through gamma-rays produced by 
annihilating H° pairs in and around the galactic cen- 
ter has been studied in several papers HH. Al- 
though IDM models (with m H o < mw and large values 
of m/j) have been shown capable of producing spectac- 
ularly clear monochromatic line signals @, the absolute 
photon fluxes are subject to large astrophysical uncer- 
tainties. Also, the most interesting photon energy range 
has not yet been covered (but is currently examined by 
the Fermi Gamma-ray Space Telescope [17[)- Thus it is 
not yet possible to derive any robust constraints from 
indirect detection. 

The WMAP data on the cosmic microwave background 
radiation limits the relic abundance of dark matter, 
Qdm, to 0.094 < QoAih 2 < 0.129 [l(|, where h is the 



Hubble constant in units of 100 km Mpc s . The lower 
bound need only be respected if the IDM alone is to solve 
the dark matter problem. However, the relic density of 
H° particles must always respect the upper bound, since 
they otherwise would contribute to too much cold dark 
matter. 

While the electroweak precision tests in the absence 
of new physics beyond the standard model favor a light 
Higgs boson, this does not remain true once the inert 
scalars are taken into account. In [J it was shown that 
the inert particles make important contributions to the 
electroweak observable T according to 



AT 



1 



24tt 2 c 



(m H ± - m A o) (m H ± - m H o) , (6) 



where a is the fine-structure constant. With appropriate 
inert scalar mass splittings the IDM can hence compen- 
sate for the otherwise too small values of T arising for 
heavy standard model Higgs bosons. From Fig. 1 in Q 
we have estimated the allowed AT range (at 68% C.L.) 
as a function of mt- Together with Eq. ((6|) this translates 
into the constraint 



(m,H± — mjsfi) {m H ± — mjjo) 



24Tr 2 av 2 



( mh 



0.15 In I — 1 ±0.1 



(7) 



As first noted in 0] and more explicitly argued in [j| 
the precise LEP I measurements of the Z boson width 



forbids the decay channel Z 
quired that 



H°A°. Hence it is re- 



(8) 



When it comes to limits set by the LEP II experiments, 
previous studies have provided rough estimates based on 
existing analyses for neutralinos, Xi> within the MSSM. 
By simply applying the same upper limits on the e + e~ — > 
H°A° cross section as those set by the different LEP 
collaborations on X1X2 production, one finds that the 
LEP II data seem capable of ruling out important parts 
of the parameter space [3, [H[ ■ These observations bring 
up the question of what restrictions a more specific study 
of H°A° production could put on the IDM parameter 
space. 



III. LEP II ANALYSIS 
A. Motivation 

Although a complete analysis including optimization of 
cuts, a detailed detector simulation and standard model 
background event generation would put the most accu- 
rate LEP II limits on the IDM, it also demands a lot of 
effort to be carried out. On the other hand, the earlier 
estimates of LEP II constraints [J, [l3[ are quite coarse 
and, as motivated below, need to be improved (or pos- 
sibly verified). In this work we have therefore chosen to 
reuse cuts, detector acceptances, simulated backgrounds 
and derived production cross section limits from an ex- 
isting MSSM analysis [l(|, and instead focus on IDM sig- 
nal event generation and efficiency determination. This 
method, which will be described in detail in Section III 
B, is careful enough to put more accurate limits on the 
H° and A masses than the current existing estimates.. 

The different LEP experiments have separately 
searched for neutralinos via the production process 
e + e~ — > y?Y ° followed by the decay x° — * Xiff 
[l8l. [T9L I20I |2lf . The lightest neutralino, Xi> is assumed 
to be stable and would show up as missing energy, while 
the fermion pair, //, could be detected as acollinear jets 
or leptons. 

As noted in [f| , there are several reasons why the lim- 
its derived for X1X2 production need not apply to H°A° 
production. The inert particles are scalars while the 
neutralinos carry spin 1/2, and in principle spin cor- 
relation effects could make a difference. The produc- 
tion of H°A° pairs only proceeds via s-channel spin-1 
Z bosons, e + er — > Z — > H°A°, forcing the spinless out- 
going scalars into states with large transverse momen- 
tum, while neutralinos have the ability to conserve an- 
gular momentum with the help of their intrinsic spin. 
Moreover, neutralino production differs in that it can 
also take place via t-channel sclcctron exchange. Simi- 
larly. A only decays according to A — > H°Z — > H°ff, 
while both xl -> x\Z -> x?// and x° 2 -> // - X°iff 
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FIG. 1: Representative Feynman diagrams contributing to 
X1X2 events at LEP. (a)-(d) show the process factorized into 
e + e~ — > X1X2 production, (a) and (b), and subsequent x° —> 



Xiff decay, (c) and (d). (e)-(h) show the unfactorized process 

xlxlfl 



are allowed (where / denotes a sfcrmion). Figs. QJa)-(d) 
and Figs. HJa)-(b) show Feynman diagrams contributing 
to the production and decay processes for the MSSM and 
IDM, respectively. In summary there are many effects 
which threaten to spoil a simple usage of the LEP II 
MSSM limits directly onto the IDM. 

To illustrate observable outcomes we take help of the 
event generator MadGraph/MadEvent [22j|, which is 
to be further discussed in Section III B. A comparison 
between the angular distributions of e + e~ — > X1X2 an d 
e + e~ — > H a A° indicates that there are significant dif- 
ferences between the two processes. Fig. [3] shows the 
angular probability distribution of produced x° arLC1 ^° 
for representative models with various neutralino and in- 
ert scalar masses. As can be seen, A is mainly produced 
with large transverse momentum, while the x° distribu- 
tion is closer to isotropic. 

Once including the decay processes one would expect 



FIG. 2: Feynman diagrams contributing to H°A° events at 
LEP. (a) and (b) show the process factorized into e + e~ — > 
H°A° production, (a), and subsequent A — » H°ff decay, 
(b). (c) shows the unfactorized process e+e" -» H°A° 
H°H°ff. 



the angular distribution of final state particles to get a bit 
smeared out. In order to properly include spin correla- 
tion and off-shell effects the full e + e~ — > X1X2 ~> XiXiff 
and e + e~ — > H°A° — ► H°H a ff matrix elements, shown 
in Figs. [IJe)-(h) and Fig. [He), have to be calculated 
without factorizing them into production (e + e _ — > XiX.2? 
e+e" -> H°A°) and decay {xl Xiff, A H°ff) 
parts (see, e.g., |23|). However, the analysis of [3 does 
not include any spin correlation effects and we hence fo- 
cus on the factorized e + e~ — > X1X21 X2 ~^ Xiff process, 
Figs. QJa)-(d), for the neutralinos. Even if spin corre- 
lations are absent for scalars, off-shell effects could in 
principle make some difference and we hence use the un- 
factorized process, Fig.^c), for the IDM. 

Fig. |4] shows the angular probability distribution of 
the final state fermions for the same models as in Fig. [3] 
Since we factorize the MSSM processes, and since the 
inert scalars carry no spin, an isotropic decay in the rest 
frame of the mother particle (\2 or ^ s expected. This 
makes the final angular distributions considerably flatter, 
and the MSSM and IDM models typically produce rather 
similar results, see Figs.gfb)-(c). However, if the velocity 
of the mother particle is large and the energy injected into 
the // pair during the decay is relatively low, the boost 
can preserve some significant difference. In other words, 
the MSSM and IDM final fcrmion distributions will still 



5 



Angular distribution of A vs x\'- mi =60, rri2=72 Angular distribution of A° vs x\'- m i=60, m^=X24 Angular distribution of ^4° vs X2 : mi =40, 7/72=160 
1 



8 °- 8 

"D 
"D 

^0.6 

C 

CD 

^0.4 



-O 

2 0.2 



IDM 


(a) 


MSSM 


s. / 
/ 


\ y 


S / 


\ ✓ 











■ IDM 
■MSSM 



(b) 




IDM 


(c) 


MSSM 


s / 


\ S 


N, / 
>» / 









-0.5 0.5 

cos 9 



-0.5 0.5 

cos 



-0.5 0.5 

cos 



FIG. 3: Angular distribution of A (IDM) versus x° (MSSM) produced in e + e — > X1X2 an d e + e — > H°A°, respectively. The 
header of each subfigure displays (in units of GeV) the mass mi of H° and Xi, an d the mass mi of A and X2- The beam pipe 
is denned to be along cosf? = ±1 and the center-of-mass energy is ^J~s~ = 206 GeV. The large difference between the IDM and 
the MSSM models is due to the scalar versus fermion nature of the outgoing states. 
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FIG. 4: Angular distribution of final state fermions from e+e" H°A° -> H°H°ff (IDM) versus e+e" X1X2, X° -* Xiff 
(MSSM). The models are the same as in Fig. [3] In (a) the velocity of the mother particle A /xi is large and the energy injected 
into the fermions during the decay is relatively low, and hence the discrepancy from Fig. [3] can survive. 
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FIG. 5: Fermion opening angle distribution for the same processes and models as in Fig. [4] /3 is the angle between the two 
outgoing fermions / and / as measured in the lab frame. In (a) the IDM and SUSY distributions happen to be overlapping, 
but one should keep in mind that in the MSSM the decay distribution can depend on the gaugino fraction of x§- Figs. 13151 are 
based on 10 5 events per model, as generated with MadGraph/MadEvent [22 |. 
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differ for models where the center-of-mass energy, y/s, 
is well above the kinematical limit for neutralino/inert 
scalar pair production, while at the same time the mass 
difference between the produced particles is small, see 
Fig.[4ja). Noteworthy, this is exactly the interesting dark 
matter preferred region of rather light H° (m#o < mw) 
together with a considerable amount of coannihilation 
with A (Am<10 GeV) $MM- 

It should be stressed that there are of course other 
observables than individual scattering angles which may 
differ between the MSSM and the IDM models. As an 
example Fig. [5] shows the // opening angle distributions, 
where clear differences are seen in Figs. [SJb)-(c). Which 
observables are important depends on what cuts are im- 
posed, and in general it can be hard to extract any reli- 
able information without doing a more complete analysis. 
In any case, it should be clear from the above examples 
that there exist differences between the MSSM and the 
IDM which could hinder a direct application of the ex- 
isting LEP II neutralino production cross section limits 
on the inert scalars. We therefore find it well motivated 
to study the LEP II bounds on the IDM more closely. 



B. Method 

1. General outline 

A DELPHI collaboration search for pair production of 
neutralinos in the center-of-mass energy range 192-208 
GeV [lH forms the basis of this paper. In [l(| the absence 
of an excess above the standard model background in the 
DELPHI data was used to put upper limits on neutralino 
production cross sections. The analysis included a care- 
ful standard model background generation and a detailed 
detector simulation. The signal events were produced 
with the help of Susygen 2.2004 and Jetset 7.4 
[25l ] was used for quark fragmentation. Every event was 
tested against a number of searches with cuts optimized 
for different final state topologies, each using a sequential 
cuts approach. These searches, which are described in de- 
tail in [16j, were ordered according to 1) acoplanar jets, 
2) acoplanar leptons, 3) multijets, 4) multilcptons, and 
5) asymmetric taus, and an event being accepted by one 
of the selections was explicitly rejected by those appear- 
ing later in the list. However, the searches were designed 
to be mutually exclusive in order to substantially reduce 
the possibility of events being capable of passing more 
than one of the searches. Based on the resulting efficien- 
cies for simulated signal events, 95% C.L. upper limits 
on neutralino production cross sections were derived and 
presented as functions of neutralino masses in contour 
plots. Of interest for this paper are Figs. 12 and 13 in 
[10 ] where X1X2 production followed by X2 decay into Xi 
plus a fermion pair (qq, [i + [i~ , e + e~ , or according to the 
branching ratios of the Z boson) was assumed. Direct 
translation of these kinds of plots is what previously has 
been used to estimate rough limits on the IDM. 



The basic strategy of this work is simple: simulate 
e + e~ — > H°A° events, pass them through the same cuts 
as those used in [l6[ to determine IDM signal efficien- 
cies, and finally rescale the upper limits on production 
cross sections derived in fl6| in accordance with the ra- 
tio between the MSSM and IDM efficiencies. However, 
in reality the full analysis of [l(| is not practically possi- 
ble to reproduce in full detail. Instead we rely on slightly 
different codes for event generation, a rough detector sim- 
ulation, and a somewhat simplified set of cuts. Thus a 
reevaluation of also the MSSM efficiencies is required, 
both for consistency and in order to check the reliability 
of our method. 



2. MSSM framework 

The free parameters of the R-parity conserving MSSM 
considered by (l6j are tan/3, fi, M2, ttlq, A and tua- The 
gaugino mass parameters Mi and M2 are assumed to be 
related according to the mSUGRA unification relation: 
Mi w O.5M2. By choosing a high common scalar mass 
parameter mo ~ 1 TeV, [16[ put focus on models with 
sfermions much heavier than the produced neutralinos. 
For fixed values of tan/3 (i.e. the ratio between the vac- 
uum expectation values of the two Higgs doublets), [l6| 
performed scans over M2 and the Higgs mass parame- 
ter fi to find models with suitable values of the x? and 
X2 masses. The common trilinear coupling A and the 
pscudoscalar Higgs mass tua typically have no noticeable 
influence on the results. 

Because of the assumed unification relation between 
Mi and M2, there are no models with the mass of 
X® much larger than twice the mass of Xi- Therefore 
1 1 61 ] considered X1X3 production followed by the decay 
X3 — > Xiff when they explored larger neutralino mass 
differences. However, their notation was to not explic- 
itly write out any X3, but instead use x° f° r anv heav- 
ier neutralino. We follow the same convention here. 
Although one can in general find several points within 
the (tan/3,M 2 ,/z) parameter space giving the same set of 
(XiiXi) masses, the differences in efficiencies are typically 
small. 



3. Signal generation 

For the generation of signal events we use the Mad- 
Graph/MadEvent package [22|. There are two main 
reasons for this choice: 1) It is possible to generate both 
MSSM and IDM events with MadGraph/MadEvent. 
2) The code allows multiparticle final states, giving the 
opportunity to generate the complete 2 — ► 4 processes in 
one step, thus fully taking spin correlation and off-shell 
effects into account. In agreement with [l6[, the MSSM 
processes are factorized into e + e~ — > X1X21 followed by 
cither xl XxQQ, X2 ~> XiH + V~, or xl Xie+e". The 
decays are arranged to be isotropic in the rest frame of 
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after which the decay products are boosted into the 
lab frame. The IDM events, on the other hand, are gen- 
erated via the unfactorized processes e + e~ — ► H°A° — > 
H {) H a ff, where // represents either qq, or e + e~. 

As a consistency check we also study numerous unfactor- 
ized MSSM processes as well as factorized IDM processes. 
Although slight deviations in for example the angular 
distribution of produced particles can be found between 
the factorized and unfactorized MSSM processes, the dis- 
crepancies left at the level of efficiency determination are 
quite small. For the IDM this agreement is even better. 

We fix the center-of-mass energy y/s to 206 GeV in our 
analysis, and note that the final results hardly depend on 
the exact value. 

Pythia 6.4 [26| is used for the fragmentation and 
hadronization of final state quarks, and jets are recon- 
structed with the PYCLUS algorithm. The minimum 
transverse jet separation is set to dj i n = 10 GeV, while 
at the same time requi ring at least two reconstructed jets, 
all in agreement with [16j |. 



4- Cuts 

Different cuts are imposed depending on whether the 
final state includes quark or lepton pairs. We mimic the 
cuts of the acoplanar jets search and the acoplanar lep- 
tons search in [lj| as closely as possible, and apply them 
to the quark and lepton events, respectively. Muon and 
electron identifications are performed with the help of 
tabulated acceptances extracted from an earlier detailed 
study of the DELPHI detector [27|. Given the four- 
momentum of the lepton, a probability for identification 
is assigned and used to generate identification status. 

Because of e.g. our lack of a detailed detector simu- 
lation, some of the cuts have to be somewhat simplified, 
and some are not even possible to carry out. However, 
much effort has been put into accounting for all major 
effects in order to assure the reliability of our final con- 
clusions. 

In principle an event not passing the acoplanar jets or 
acoplanar leptons cuts could instead happen to be ac- 
cepted by one of the subsequent searches (i.e. the multi- 
jets, multilcptons, or asymmetric taus selections). Owing 
to the design of the searches we expect these events to 
be quite rare, and that the relative MSSM and IDM ef- 
ficiencies in any case can be fairly represented by the 
corresponding results from the acoplanar searches. We 
therefore only impose one set of cuts on each event: the 
acoplanar jets (for qq final states) or the acoplanar lep- 
tons (for fi + fi~ and e + e~ final states) selections. The 
exact selection cuts are described in appendices A and 
B. 



5. Calculating efficiencies 

We start our analysis by generating (on average a 
handful of different) MSSM models for each set of (X11X2) 
masses considered. We create the models using the spec- 
trum generator Suspect 2.34 [HI, and feed them into 
MadGraph/MadEvent to simulate MSSM events. In 
total we generate more than 15 x 10 6 X1X2 events dis- 
tributed among roughly 200 different MSSM models. Ef- 
ficiencies are then calculated by passing the generated 
events through our cuts. As a reliability check of our 
analysis, the found efficiencies are compared with those 
of Table 8 in [lj|- Although an exact agreement is too 
much to demand, the ratios between our calculated ef- 
ficiencies and those in [l6[ should at least stay rather 
constant for all (xiiX.2) niass points. Indeed we find that 
this ratio stays in the interval 0.7-1.0 for most models, 
a result accurate enough for our needs. Some larger de- 
viations occur, but this happens mainly in not very im- 
portant regions of the parameter space, and can possibly 
be understood by taking into account the uncertainties 
in the jet energy determination by the true detector. See 
Section V for further discussion. 

Within the IDM the cross sections for production of 
H°A° events are on the other hand completely speci- 
fied once rune and m^a are given, and hence only one 
model for each set of (H°,A°) masses needs to be exam- 
ined. Again MadGraph/MadEvent is used to simu- 
late (about 4 x 10 6 ) signal events, which are subsequently 
taken through the same cuts in order to extract the cor- 
responding IDM efficiencies. 



6. Constraining the IDM 

We can now determine the ratio between our calculated 
MSSM and IDM efficiencies as a function of the produced 
ncutralino/inert scalar masses. These are then used to 
rescale the cross section upper limits presented in [l6j]. 
It should be noted that although the 14 mass combina- 
tions given in Table 8 in [l6j serve as our starting check 
points, we investigate many more to cover the complete 
(mjjo ,m^o) plane and to be more detailed around more 
critical regions. 

In the decay of ^4°, fermion pairs are produced at tree 
level via s-channel Z bosons, and the branching ratios 
into fermions are therefore expected to be close to those 
of an on-shell Z boson. It is therefore natural to fo- 
cus on the cross section limits provided by Fig. 13(d) 
in [X 61 ] . where the final fermions arc assumed to be pro- 
duced in accordance with the branching ratios of the Z 
boson. This is clearly appropriate for models with a large 
mass splitting Am (and specifically when Am > mz so 
that Z can be on-shell), but some modifications may be 
needed for models where the intermediate Z boson is 
very virtual and the branching may be somewhat dif- 
ferent. In other words, when Am is small we need to 
evaluate the branching ratios to find out whether some 
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FIG. 6: Production cross section upper limits as extracted 
from Fig. 13(d) in [l(|. For models inside the (red) solid 
[(green) dashed] contour the limits are rescaled by a factor 
0.9 (1.1) before being applied to H°A° production. The solid 
(dark blue) contour lines indicate the e + e~ — ► H°A° cross 
section. The (red) dotted-shaded region, where m H o+ m A o < 
mz, is excluded by LEP I data on the Z boson width. The 
upper right dashed line shows the LEP II kinematical limit. 



modifications might be needed (this calculation is done 
with MadGraph/MadEvent). 

Finally, we calculate the H°A° production cross sec- 
tion as a function of m#o and m A o , and compare it with 
our derived cross section upper limits in order to con- 
strain the IDM parameter space. 

IV. RESULTS 

Under our imposed cuts the resulting IDM and MSSM 
efficiencies turn out to be quite similar, an appealing, 
although not at all trivial, result. 

The efficiencies arc first determined for each individ- 
ual channel (qq, e + e _ ), after which those are com- 
bined into an efficiency representing the actual branching 
ratio. This combination is done by weighting the chan- 
nels in accordance with the decay branching ratios of 
the Z boson (i.e. the qq efficiency is given the highest 
weight). 

In general we observe that the ratio between our de- 
rived IDM and MSSM efficiencies is quite insensitive to 
the very details of the imposed cuts, and we estimate our 
sensitivity in determining this ratio to be of the order of 
10%. 

We find that whenever m H o < 80 GeV the IDM effi- 
ciencies typically are a few percent higher than those of 
the corresponding MSSM models. An important obser- 
vation is that we find no mass combinations in this region 
where the MSSM gives a higher efficiency than the IDM, 
and it is therefore appropriate to apply at least as hard 
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FIG. 7: LEP exclusion plot. The (red) dotted-shaded region 
indicates the region of the (m H o ,m A o) plane excluded by LEP 
data. The lower left triangle, where m H o + m A o < mz, is 
excluded by LEP I data on the Z boson width. The remaining 
part of the shaded region is excluded by our LEP II analysis. 
Shown is also the LEP II kinematical limit. Since we are 
assuming m H o < m A o the upper left region is not accessible. 

production cross section upper limits on the inert scalars 
as those put on the neutralinos in [l6[. 

In the specific region defined by 8 GcV< Am < 15 GeV 
and rriHO < 85 GeV, the IDM efficiencies are found to be 
about a factor 1.15-1.20 higher than those of the MSSM. 
On noting that the models with the lowest Am have a 
slightly higher branching into neutrinos compared to or- 
dinary Z boson decay, we in this region adopt a conser- 
vative factor of 0.9 with which we rescale the neutralino 
production limits given in Fig. 13(d) in [l6| . This region 
is encircled with a (green) dashed line in Fig. [5J 

Among the remaining m^o > 80 GeV models we find 
some for which the ratio between the IDM and MSSM 
efficiencies drops down to 0.9. We therefore use a factor 
of 1.1 for the rescaling here, and this region is encircled 
with a red solid line in Fig. [6j 

Except for in the low Am and high m#o regions men- 
tioned above we find it appropriate to apply the same 
production limits as for the neutralinos. While this might 
be argued to be too conservative, the points where harder 
limits could possibly be imposed are anyway far from ex- 
cluding any IDM model. 

By utilizin g t he limits on the X1X2 production from 
Fig. 13(d) in [16| we find, after rescaling, upper limits on 
the H°A° production cross section as a function of m#o 
and m A o. The cross section limits, and the regions where 
we impose rescaling, are found in Fig. [5] Comparing 
these with the calculated e + e~ — > H°A° cross sections, 
which also are shown in Fig. [6l finally tells us which IDM 
models are excluded. 

The resulting exclusion plot is shown in Fig. [7] 
Roughly speaking, our LEP II analysis exclude models 
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FIG. 8: LEP II limits and dark matter, m h = 200 GeV. The 
LEP excluded region of Fig. [7J together with a (green) gray 
region indicating points in the (m H o ,m A o) plane where there 
exist models with mh = 200 GeV capable of providing a good 
dark matter (DM) candidate H° with relic density in agree- 
ment with WMAP. Models above the (green) gray region are 
still allowed but cannot account for all dark matter. Models 
below the (green) gray region leave a relic density higher than 
the WMAP upper limit and are hence ruled out. 



FIG. 9: LEP II limits and dark matter, m h = 500 GeV. The 
LEP excluded region of Fig. [7J together with a (green) gray 
band indicating points in the (m H o ,m A a) plane where there 
exist models with mh = 500 GeV capable of providing a good 
dark matter (DM) candidate H° with relic density in agree- 
ment with WMAP. Models above the (green) gray band are 
still allowed but cannot account for all dark matter. Models 
below the (green) gray band all leave a relic density higher 
than the WMAP upper limit and are hence ruled out. 



satisfying m#o < 80 GeV, m^o < 100 GeV and Am > 8 
GeV. The sharp transition at A?ti = 8 GeV comes from 
the steep gradient of the cross section upper limit present 
in Fig. 13(d) in 0. 

In [J] it has previously been estimated that the H°A° 
production cross section is below the existing LEP II up- 
per limits for models with m H a ss 70 GeV and Am < 10 
GeV. We note that our results put somewhat harder 
constraints in that region. Moreover, fl3T ] presented 
m H o + m A o > 130 GeV as a rule of thumb for mod- 
els to be consistent with the LEP II data. For m H o > 30 
GeV our analysis hence imposes harder constraints than 
that. 

In order to find out what consequences this LEP II 
constraint has for H° as a dark matter candidate we per- 
form an extensive scan over the IDM parameter s pac e 
using DarkSusy [29[ interfaced with FormCalc [30T|. 
The relic density is calculated, including coannihilations 
with A and H ± , for models respecting Eqs. |g|), 
and ((HI). We also include the limit m H ±> 70 GeV, which 
is based on direct translation of the limits from LEP II 
chargino search results [l3| . 

Fig. Hand Fig. 1 show, for m /l = 200 GeV and m h = 500 
GeV, respectively, the regions where H° can have the 
right relic abundance to account for all dark matter. 

Because the only existing (tree level) H° self- 
annihilation channel into fcrmions is via s-channcl h ex- 



change, models with a heavy h can only provide a relic 
density in agreement with the WMAP data if either: 1) 
Am is small so that coannihilation with A becomes im- 
portant, or 2) m H o ~ 75 GeV so that the efficient an- 
nihilations into massive gauge bosons are of just about 
the right magnitude in the early Universe. If the coanni- 
hilations are too strong, or if m#o is too high, the relic 
density will fall below the WMAP lower limit, in which 
case the model is still allowed but cannot explain all dark 
matter. These regions are easily identified in Fig. [9l 

In Fig. [5] the mass of h is small enough to allow the H° 
self-annihilating cross section into fermions to be large 
enough to obtain a correct relic density in a larger part 
of the (m H o ,m^o) plane. The smaller value of mh also al- 
lows, because of cancellations between the diagrams con- 
tributing to annihilations into massive gauge bosons, for 
a slightly heavier dark matter candidate H°. 

As can be seen in Figs. [5] and Fig. [5] the LEP II data 
can exclude many IDM models of interest from a dark 
matter perspective, although definitely not all of them. 

V. DISCUSSION 

Although we find that appropriate upper limits to im- 
pose on the production cross section for inert scalars are 
similar to those already derived for neutralinos in [l6|, 
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this result is a priori far from trivial owing to the differ- 
ent properties of the particles. 

It is encouraging that the effect of including realistic 
selection efficiencies is not to invalidate, but rather to in- 
crease the confidence in simple applications of limits from 
other MSSM searches directly onto the IDM. However, 
one must keep in mind that rather different sets of cuts 
may be used in other searches, and hence any robust pre- 
diction requires an evaluation of the corresponding IDM 
efficiencies. 

One may worry that some of the simplifications used 
in this analysis misses some important effects. For ex- 
ample we do not include any initial state radiation from 
the incoming e + e~ pair. While this is certainly impor- 
tant for the models excluded by LEP I which satisfy 
mjjo + m^a < mz (and for which the intermediate s- 
channcl Z bosons can be put on-shell by radiating off 
photons from the incoming particles) we do not expect 
the initial state radiation to be very important for our 
mass range of interest. 

Another concern is that the detector generally fails to 
detect some amount of the energy in high energy jets. 
On assuming a perfect jet energy determination we in- 
deed find that for models with Am > 100 GeV our derived 
MSSM efficiencies do not match those of [lj| as well as 
they do for lower Am models. However, by artificially 
scaling down the jet energy of each event with on av- 
erage 10%-20% before applying the cuts we find good 
agreements also for these models. An important obser- 
vation to keep in mind is that we find the ratios between 
the IDM and MSSM efficiencies to be practically insen- 
sitive to the downscaling magnitude. Therefore our final 
IDM limits do not depend on the very details of the true 
energy losses. Also, the H°A° production cross section 
is in general too low within the high Am region for any 
models to be excluded anyway. 

In [l6| the acoplanar searches were each divided into 
four subselections (these subselections are the last l)-4) 
cuts presented in each of appendices A and B). Within 
their analysis these were treated as independent selec- 
tions, and we should hence actually check the efficiencies 
under each of them separately when trying to extract 
the IDM limits. However, the ratios between our derived 
IDM and MSSM subsclcction efficiencies are found to be 
practically equal to the values for the corresponding com- 
plete acoplanar selections. We hence conclude that this 
does not constitute a problem for our method. 

While our analysis certainly includes some approxima- 
tions, we have consistently been quite conservative in or- 
der to minimize the risk of overestimating the final con- 
straints on the IDM. Hence, our derived LEP II limits on 
the inert scalar masses should be robust. 



VI. SUMMARY 

In this paper we have investigated what limits can 
be inferred on the IDM from an already existing LEP 



II neutralino search. In performing this translation we 
have used a method respecting the differences between 
inert scalars and neutralinos. We have shown that the 
LEP II data can exclude a significant part of the param- 
eter space, but also that many models providing a good 
WIMP dark matter candidate are still valid. 
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Appendix A 

Below follows the acoplanar jets selection cuts we im- 
posed on the quark final states after hadronization in 
Pythia and subsequent rejection of particles lost in the 
beam pipe (which is taken to cover a polar angle of 10°). 
A particle was assumed to carry time projection cham- 
ber (TPC) information if and only if it was more than 
25° away from the beam axis. (See [3l[ for information 
about the TPC pad rows in the DELPHI detector.) Al- 
though some simplifications have been unavoidable, most 
of our cuts are very close to the original ones. For a list 
of original cuts see [l|| . 

• A minimum of two charged particles, with at least 
one of them having a transverse momentum above 
1.5 GeV, and a total transverse energy (defined as 
the sum of the absolute values of the transverse 
momenta of all individual particles) above 4 GeV, 
was required. 

• At least five charged particles had to carry TPC 
information. 

• The scalar sum of momenta of the particles carrying 
TPC information had to exceed both 4 GeV and 
10% of the total jet energy. 

• Exactly two jets were required. 

• All jets had to have a polar angle above 10°. 

• Each jet had to contain at least one particle with 
TPC information. 

• The particles emitted within 30° of the beam axis 
had to carry less than 60% of the total jet energy. 
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The absolute value of the cosine of the polar angle 
of the total momentum had to be smaller than 0.9. 

The total transverse momentum (defined as the 
norm of the vector sum of the individual transverse 
momenta) had to be larger than 6 GeV. 

The momentum of a jet divided by its energy had 
to exceed 0.5 for the most energetic jet and 0.4 for 
the second most energetic jet. 

The total transverse momentum had to be larger 
than 15 GeV, or the acollinearity (defined as the 
supplement of the angle between the jets) had to 
exceed 40°, or the ratio between the momentum 
and the energy of the second most energetic jet had 
to be above 0.8. 

If the total transverse momentum was below 12 
GeV, the average momentum of the particles with 
TPC information had to lie between 0.8 GeV and 8 
GeV, and the energy of the most energetic neutral 
particle had to be below 40% of the total jet energy 
(if the total jet energy exceeded 20 GeV) or below 
5 GeV (if the total jet energy was lower than 20 
GeV). 

If the total jet energy was below 50 GeV, the aver- 
age momentum of the particles with TPC informa- 
tion had to lie between 0.8 GeV and 8 GeV, and 
the momentum of the most energetic charged par- 
ticles had to lie between 5% and 70% of the total 
jet energy. 

If the total jet energy was below 20 GeV, the aver- 
age momentum of the particles with TPC informa- 
tion had to lie between 0.8 GeV and 4 GeV but at 
least be larger than 20% of the total jet energy, the 
momentum of the most energetic charged particle 
had to lie between 10% and 60% of the total jet 
energy, and the energy of the most energetic neu- 
tral particle had to be below 35% of the total jet 
energy. 

If there were any neutral particles, the norm of the 
vector sum of the jet momenta and the momentum 
of the most energetic neutral particle had to exceed 
2.5 GeV. 

All neutral particles had to have energies less than 
60 GeV. 

All charged particles had to have energies less than 
20 GeV. 

The jets had to pass one of the following cuts: 

1. The invariant mass had to be below 10% of 
y/s, the missing mass had to be above 70% 
of y/s, the total transverse momentum had to 
be above 7 GeV, and the scaled acoplanarity 



(defined as the sine of the minimum angle be- 
tween a jet and the beam axis times the sup- 
plement of the angle between the projections 
of the jets onto a plane perpendicular to the 
beam axis) had to exceed 40°. 

2. The invariant mass had to lie between 10% 
and 30% of y/s, the missing mass had to be 
above 60% of y/s, the total transverse momen- 
tum had to be above 8 GeV, and the scaled 
acoplanarity had to exceed 25°. 

3. The invariant mass had to lie between 30% 
and 50% of y/s, the missing mass had to be 
above 45% of y/s, the total transverse mo- 
mentum had to lie between 12 GeV and 35 
GeV, the total longitudinal momentum (de- 
fined as the norm of the vector sum of the 
individual longitudinal momenta) had to be 
below 35 GeV, the scaled acoplanarity had to 
exceed 25°, and the acollinearity had to be 
lower than 55°. 

4. The invariant mass had to lie between 50% 
and 70% of y/s, the missing mass had to be 
above 20% of y/s, the total transverse momen- 
tum had to lie between 12 GeV and 35 GeV, 
the total longitudinal momentum had to be 
below 35 GeV, the scaled acoplanarity had to 
exceed 25°, and the acollinearity had to be 
lower than 55°. 



Appendix B 

Below follows the acoplanar leptons selection cuts we 
imposed on the lepton final states. A particle was as- 
sumed to carry TPC information if and only if it was 
more than 25° away from the beam axis. (See [3l[ for 
information about the TPC pad rows in the DELPHI 
detector.) Although some simplifications have been un- 
avoidable, most of our cuts are very close to the original 
ones. For a list of original cuts see [l6j . 

• Two leptons needed to be identified. 

• At least one of the leptons had to have a transverse 
momentum above 1.5 GeV, and a total transverse 
energy above 4 GeV was required. 

• Each lepton had to have an energy above 1 GeV. 

• Each lepton had to carry TPC information. 

• The acoplanarity between the leptons had to be 
larger than 10°. 

• The acollinearity between the leptons had to be 
larger than 10°. 

• The absolute value of the cosine of the polar angle 
of the total momentum had to be smaller than 0.9. 
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• The total transverse momentum had to be larger 
than 6 GeV. 

• The leptons emitted within 30° of the beam axis 
had to carry less than 70% of the total lepton en- 
ergy. 

• If the total transverse energy was below 100 GeV 
and the total missing momentum was above 45 
GeV, then the charge times the cosine of the polar 
angle of the most energetic lepton was required to 
be positive. 

• If the momentum of the most energetic lepton was 
between 40 GeV and 100 GeV and the total energy 
within 15° of this lepton was between 40 GeV and 
100 GeV, then the charge times the cosine of the 
polar angle was required to be greater than -0.65 
for both leptons. 

• The leptons had to pass one of the following cuts: 

1. The invariant mass had to be below 10% of 
y/s, the missing mass had to be above 70% of 
y/s, the total transverse momentum had to be 



above 7 GeV, and the acoplanarity (defined 
as the supplement of the angle between the 
projections of the leptons onto a plane per- 
pendicular to the beam axis) had to exceed 
40°. 

2. The invariant mass had to lie between 10% 
and 30% of y/s, the missing mass had to be 
above 45% of y/s, the total transverse momen- 
tum had to be above 10 GeV, and the acopla- 
narity had to exceed 25°. 

3. The invariant mass had to lie between 30% 
and 55% of y/s, the missing mass had to be 
above 20% of y/s, the total transverse momen- 
tum had be above 12 GeV, and the acopla- 
narity had to exceed 15°. 

4. The invariant mass had to lie between 55% 
and 70% of y/s, the missing mass had to be 
above 20% of y/s, the total transverse momen- 
tum had to be above 12 GeV, and the acopla- 
narity had to exceed 15°. 
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